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We study the hot electroweak phase transition (EWPT) by 4-dimensional lattice 
simulations on lattices with symmetric and asymmetric lattice spacings and give 
the phase diagram. A continuum extrapolation is done. We find first order phase 
transition for Higgs-boson masses mjj < 66.5 ± 1.4 GeV. Above this end point a 
rapid cross-over occurs. Our result agrees with that of the dimensional reduction 
approach. It also indicates that the fermionic sector of the Standard Model (SM) 
may be included perturbatively. We get for the SM end point 72.4 ± 1.7 GeV. 
Thus, the LEP Higgs-boson mass lower bound excludes any EWPT in the SM. 



1 Introduction 

The observed baryon asymmetry is finally determined at the EWPT 0. To 
understand this asymmetry a quantitative description of EWPT is needed-. 
Since the perturbative approach breaks down for large Higgs-boson masses □ 
(e.g. rriH > 70 GeV) lattice MC simulations are necessary. (Another, partly 
analytic approach has been presented in El.) 

Previous works show that the strength of the EWPT gets weaker as the 
Higgs-boson mass increases. The line of the first order phase transitions, sep- 
arating the symmetric and broken phases on the m# — T c plane, has an en< 
point, rriH.c- 3D results show that for. run > 95 GeV no EWPT exists 



moreover m,H, c ~ 67 GeVB'Qu. In 4Dl3 m# iC « 80GeV was obtained. Also 
in 4D at ran ~ 80 GeV the EWPT turned out to be extremely weak, even 
consistent with the no phase transition scenario on the 1.5-cr levella. However, 
a discrepancy appeared: the 4D estimate for the end point was higher than 
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Figure 1: Imaginary part ot-^rst Lee- Yang zero as a function of A from simulations on 
symmetric lattices with Lt = 2 £3. Filled symbols are obtained without A-reweighting, while 
open symbols with A-reweighting from the filled symbol with same shape. 

the 3D estimates. This discrepancy has been resolved in00, to be reviewed 
here. 



2 End point analysis 

We study the 4D SU(2)-Higgs model on both symmetric and asymmetric 
liil lattices, i.e. lattices with equal or different spacings in temporal (at) and 
spatial (a s ) directions. In the asymmetric case equal lattice spacings are used 
in the 3 spatial directions (a^ = a s , i = 1,2, 3). The asymmetry of the lattice 
spacings is given by the asymmetry factor £ = a s /at- The different lattice 
spacings can be ensured by different coupling strengths in the action iqr time- 
like and space-like directions. The action reads in standard notationE 
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We introduce and P 2 = (3 s (3f The anisotropies 7^ = /3t//3 s and 

7k = K t/ K s are functions of £. These functions have been determined pertur- 
batively and non-perturbativelyll3 demanding the restoration of the rotational 
symmetry in different channels. We use £ = 4.052, which corresponds tp-r/K = 4 
and 7/3 = 3.919. Details of the simulation techniques can be found inl!3. 
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Figure 2: Schematic phase diagram. The solid line represents the LCP defined by the end 
point condition. The numbers on the line correspond to the temporal extension (the dashed 
lines show their projection to the re - A plane). The dotted lines running into these numbered 
points correspond to first order phase transitions for = const, but different Rhw-s. A 
LCP defined by a constant Rjjw value is shown by the long dashed line. 



The determination of the end poinL-pf the finite temperature EWPT is 
done by the use of the Lee- Yang zeros E-3 of the partition function Z. Near 
the first order phase transition point the partition function reads Z = Z s + 
Zi, oc exp(— Vf s ) + cxp(— Vfb) , where the indices s(b) refer to the symmetric 
(broken) phase and / stands for the free-energy densities. We also have /& = 
f s + a(n — k c ) , since the free-energy density is continuous, ft follows that 
Z oc exp[— V(f s + /&)/2] cosh[— Va(n — k c )] which shows that for complex k 
Z vanishes at Im(/i) = 2tt ■ (n — 1/2) /(Va) for integer n. In case a first order 
phase transition is present, these Lee- Yang zeros move to the real axis as the 
volume goes to infinity. In case a phase transition is absent the Lee- Yang 
zeros stay away from the real k axis. Thus the way the Lee- Yang zeros move 
in this limit is a_good indicator for the presence or absence of a first order 
phase transition t3. Denoting kq the lowest zero of Z, i.e. the position of 
the zero closest to the real axis, one expects in the vicinity of the end point 
the scaling law Im(/to) = G{L tl XjV~" + n%{L t , A). In order to pin down the 
end point we are looking for a A value for which Kq vanishes, pip practice we 
analytically continue Z to complex values of k by reweightingE£l the available 
data. Also small changes in A have been taken into accqimt by reweighting. 
The dependence of k§ on A from our symmetric simulations!!!] is shown in Fig. 1. 
To determine the critical value of A i.e. the largest value, where Kq = 0, we 
have performed fits linear in A to the non-negative Kq values. 

In the isotropic lattice simulation 0, we have used L t = 2. The Lee- 
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Figure 3: Dependence of Rhw,c i.e. Rhw corresponding to the end point of first order 
phase transitions, on and extrapolation to the continuum limit. 

Yang analysis gave A c = 0.00116(16) for the end point. Performing T = 
simulations with the same parameters this can be converted to a Higgs boson 
mass. The value m,H, c = 73.3±6.4GeV has been obtained, _azhich is compatible 
with our estimate based on a study of Binder cumulantsE3 and the previous 
4D studies i'l. 

In the anisotropic lattice simulation caselld we also performed a continuum 
extrapolation, moving along the lines of constant physics (LCP). A schematic 
illustration is shown in Fig. 2. 

The technical implementation of the LCP idea has been done as follows. 
By fixing (3 — 8.0 in the simulations, we have observed that g^ is essentially 
constant within our errors. For the small differences in gn we have performed 
perturbative corrections. We have carried out T ^ simulations on L t = 
2, 3, 4, 5 lattices (for the finite temperature case one uses L t <c L x , L y , L z ), and 
tuned k to the transition point. This condition fixes the lattice spacings: a t = 
a s /£ = l/(T c L t ) in terms of the transition temperature T c in physical units. 
The third parameter A, finally specifying the physical Higgs mass in lattice 
units, has been chosen (using the Lee- Yang analysis) so that the transition 
corresponds to the end point of the first order phase transition subspace. 

Having determined the end point X c (L t ) for each L t we calculated the 
T = quantities (RhWi9r) on V — (32L t ) • (8L t ) • (6L t ) 3 lattices. Having 
established the correspondence between X c (L t ) and Rhw,c the L t dependence 
of the critical Rhw,c is easily obtained. Fig. 3 shows the dependence of the 




4 



2.4 


1 






j\ symmetric- 


phase 


2.2 
a: 

B 






Higgs-phaseJ \ 




y 
E- 






2 
1.8 


~ 


\endpoint 

~ 



_J I I I I I I I I I I I I 

0.2 0.4 0.6 0.8 

R HW 



Figure 4: Phase diagram of the SU(2)-Higgs model in the (T c /m,H — Rhw) plane. The 
continuous line — representing the phase-boundary - is a quadratic fit to the data points. 

end point Rhw values on \fL\. A linear extrapolation in yields the 

continuum limit value Rhw.c = 0.83 ± 0.02, which corresponds to m,H,c = 
66.5 ±1.4 GeV. Note that m#_._decreases for increasing L tn| This observation 
resolves the discrepancy of 3Dl3ol3 and previous L t — 24DHa results. 

Comparing our result to those of the 3D analyses ErtiQ one observes com- 
plete agreement. Since the error bars on the end point determinations are on 
the few percent level, the uncertainty of the dimensional reduction procedure 
is also in this range. This indicates that the analogous perturbative inclusion 
of the fermionic sector results als£LiiLfew percent error an rriH,c- 

Based on our published dataE3'E3 and the results ofliil we draw the precise 
phase diagram in the (T c /niH — Rhw) plane of the SU(2)-Higgs model. This 
is shown in Fig. 4. The continuous line - representing the phase-boundary - 
is a quadratic fit to the data points. 

Finally, we determined the end point value in the full SM. We use per- 
turbation theory Il3 to transform the SU(2)-Higgs model end point value to 
the full SM. We obtain 72.4 ± 1.7 GeV. The dominant error comes from the 
uncertainty on the position of the end point. 

3 Conclusions 

In conclusion, we have determined the end point of hot EWPT with the tech- 
nique of Lee- Yang zeros from simulations in 4D SU(2)-Higgs model. The phase 
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diagram has been also presented. The transition is first order for Higgs masses 
less than 66.5 ± 1.4 GeV, while for larger Higgs masses only a rapid cross-over 
is expected. Our results show that integrating out the heavy modes pertur- 
batively is sufficiently precise. Thus the above value can be perturbatively 
transformed to the full SM, yielding 72.4 ±1.7 GeV for the end point Higgs 

mass. . . 

The experimental lower limit of the SM Higgs-boson mass is 89.8 GeVliH. 
Taking into account all errors, our end point value excludes the possibility of 
any EWPT in the SM. Thus outwork emphasizes the need of EWPT analyses 
based on extensions of the SMLL3. 
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